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The experimental achievements during the past year in demonstrating the existence of long-ranged spin-triplet 
supercurrents in ferromagnets proximity coupled to singlet superconductors open up the possibility for new in- 
teresting physics and applications [for a review, see M. Eschrig, Phys. Today 64(1), 43 (2011)]. Our group 
reported the injection of triplet supercurrents into a magnetically uniform ferromagnet (Co) by sandwiching 
it between two helimagnet/superconductor (Ho/Nb) bilayers to form a Nb/Ho/Co/Ho/Nb-type Josephson de- 
vice. In the function of the Ho layer thicknesses, the supercurrent was found to modulate in a complex way 
that seemed to depend on the magnetic structure of Ho. To understand this unusual behavior, we have theoreti- 
cally studied the properties of an ideal Josephson device with a helimagnet/ferromagnet/helimagnet (HM/F/HM) 
barrier in the clean limit using the Eilenberger equation; we show, in particular, that the maximum triplet super- 
current that can pass across the barrier will depend non-monotonically on the thicknesses of the HM layers if 
the HM and F layers are magnetically exchange coupled at their interface. 



I. INTRODUCTION 

The intense interest in understanding the interplay be- 
tween superconductors (S) and ferromagnets (F) was primar- 
ily triggered by the pioneering work of Ryazanov et al.} who 
showed that the supercurrent in a S/F/S junction containing 
the weak ferromagnet CuNi was periodically modulated: first 
by temperature as a result of varying the exchange energy of 
the low Curie temperature CuNi, and also in the function of 
the CuNi thickness. 2 The oscillating and decay lengths of the 
thickness modulation are much shorter (£ w 1 nm) if a strong 
ferromagnet is used, such as Co, 3-5 Fe, 6 ' 7 or Ni, 3 ' 8 but because 
of their long mean free paths for both scattering and spin flip 
they are attractive materials to use in S/F/S junctions. For a 
detailed review on the singlet proximity coupling of S and F 
materials see Ref. 9 and references therein. 

The decay lengths of supercurrents in ferromagnets can be 
radically extended if the electron pairs transform from a sin- 
glet to a triplet spin state at a S/F interface via a spin-mixing 
process. This is known as the long-range triplet proximity ef- 
fect; for a comprehensive review, see Ref. 10 and references 
therein. One way to theoretically promote the spin-mixing 
conversion between singlet pairs and triplet pairs at a S/F in- 
terface is to incorporate a magnetically inhomogenous layer 
between the S and F layers. 11 The triplet electron pairs which 
form are expected to be much more robust in ferromagnets, 
with longer coherence lengths (£ 3> 1 nm) than pairs in a sin- 
glet state. Therefore the demonstration of a supercurrent in a 
S/F/S junction where the F layer thickness exceeds any length 
scale possible for a singlet pair to exist is considered indirect 
proof of the long-range triplet proximity effect. Signatures of 
such a long-range proximity effect were first reported more 
than two decades ago (see, e.g. Ref. 12). However, it was 
only in 2006 when the first major breakthrough was achieved 
with supercurrents reported in planar Josephson devices with 
half-metallic (i.e., fully spin polarized) barriers that were hun- 
dreds of nanometers long; 13 see also the related articles in Ref. 
14. In the same year, a triplet superconducting state was also 
reported in Ho, a helimagnetic rare-earth metal, which formed 



the junction of a superconducting interferometer. 15 These re- 
sults stimulated intense theoretical work aimed at understand- 
ing triplet supercurrents better and how they could be created 
at S/F interfaces in a more routine way. 

Over the past year it would seem that the puzzle to control 
triplet pair generation in S/F/S devices was finally solved with 
long-ranged supercurrents consistent with spin-triplet theory 
demonstrated in a wide range of ferromagnets, 16 including 
a Co/Ru/Co synthetic antiferromagnet interfaced by normal 
(i.e., non-magnetic) metal spacers and ferromagnetic alloys 17 
(see also Ref. 18); a Ho/Co/Ho composite barrier 19 (see also 
Refs. 20-22); a magnetic Cu 2 MnAl Heusler compound with a 
complex magnetic profile 23 (see also Ref. 24); a half-metallic 
Cr02 wire; 25 and a Co nanowire. 26 In addition to these ex- 
periments, superconducting gap features have been very re- 
cently measured by scanning tunneling spectroscopy 27 in the 
half-metallic manganite La 2 /3Ca 1 / 3 Mn03 (LCMO) grown on 
YBa2Cu307_a. The gap features were observed in the LCMO 
up to a thickness of approximately 30 nm, and so offer the 
potential to explore the fundamental properties of a triplet su- 
perconducting state. 

The experiments reported in Refs. 17, 19, and 23 share 
many similarities with the theoretical triplet junction proposed 
by Houzet and Buzdin: 28 a S/F'/F/F"/S junction in which the 
magnetizations in the F' and F" layers are non-collinear with 
the magnetization in the F layer, allowing control over the 
creation of triplet electron pairs. Our group reported Joseph- 
son devices in which the F' and F" layers were substituted 
by Ho and the F layer used was Co. 19 The maximum super- 
current (i.e., the Josephson critical current Ic) in these de- 
vices was found to depend non-monotonically on the thick- 
ness dh of each Ho layer, with peaks in Ic when dh ~ 4.5 
nm and dh ss 10 nm. These thicknesses appeared to corre- 
spond to an optimum spin-triplet proximity effect, and by in- 
creasing the Co thickness df in the 0—16 nm range, a slow 
decay in Ic was observed with a decay length of ~ 10 nm at 
4.2 K. This decay length is almost ten times larger than that 
found in simple (spin-singlet dominated) Co-based Josephson 
devices. 3 ' 4 Alidoust and Linder 20 established theoretically that 



2 



spin-triplet supercurrents could account for the slow decay 
in Ic with df, but their model did not account for the com- 
plex dependence of Ic on d%. When rare-earth ferromagnets 
are grown on transition-metal ferromagnets, it is well known 
that the two materials magnetically couple at the interface. In 
view of this fact, we solve the Eilenberger equation for an 
ideal S/HM/F/HM/S Josephson device with arbitrary interfa- 
cial magnetic coupling at the HM/F interfaces, and demon- 
strate a profound dependence of the triplet supercurrent am- 
plitude on the thicknesses of the HM layers. The results pro- 
vide a new level of understanding into the conversion process 
between singlet and triplet Cooper pairing at a superconduc- 
tor/helimagnet/ferromagnet interface. 

II. GENERAL FORMALISM 

The S/HM/F/HM/S device considered in this paper is il- 
lustrated in Fig. 1. All layers are normal to the z direction. 
The HM layers have a thickness dh and their magnetization 
rotates in the {x, y} plane as a function of z? 1 The strongly 
ferromagnetic F layer has a thickness d t and a uniform mag- 
netization in the +x direction. The S layers at z < and 
z > d = 2dh + df are infinitely thick and in a singlet state. 
Subscripts h, /, and s refer to the HM, F, and S layers. 

We assume a moderately clean limit and a temperature 
close to the critical temperature Tc of the superconducting 
leads, therefore we adapt the linearized Eilenberger equa- 
tion (LEE). Due to the presence of non-uniform magnetiza- 
tion, we consider both singlet and triplet correlations. The 
anomalous Green's function is then a 2 x 2 matrix given by 
l(z,6,u) = fl + f-l where / = (/- , p ',/*), and 
<t = (& x , a y ,a z ) is a vector of the Pauli matrices. The com- 
ponent /° is associated with spin-singlet pairs, while the three 
components in / are associated with spin-triplet pairs. To 
understand the meaning of the different triplet components, 
we switch to an alternative representation where we have 
/ = f°i + f x a x + f + &+ + f-&- withfj* = (&y±ia z )/V2 
and /± = (p =p if z )/V2- The components f x and /± then 
correspond to triplet pairs with spin projections and ±1 to 
the x direction, respectively. 

Since the HM layers have a weak magnetization, all compo- 
nents in / are non-zero. The Fermi surface average is (/) « 
in the moderately clean limit 29 and so the LEE reads 

v h cos 9 ^ + (2w + r^ 1 ) h + il h ■ 0, f h } = 0, (1) 

where {x, y} = xy + yx is the anticommutator, while v a 
and T a are the Fermi velocity and the pair-breaking time 
in a generic layer a. The electron mean free path is then 
@a = v aT a - The Matsubara frequencies are given by to — 
7r/csT(l + 2n)/h, and 1^ is the magnetic exchange field 
in units of frequency. This field takes the form Ih — 
Ih{cos[Qz + a],sin[Qz + a),0) in the HM layers, where 
z = \z — d/2\ — df/2is the distance to the nearest HM/F inter- 
face and Q is the wave vector of the magnetic helix. The mag- 
netic exchange coupling at the HM/F interfaces is included by 




FIG. 1: (Color online) (A) Illustration of a superconducting 
losephson device with a helimagnet/ferromagnet/helimagnet barrier 
(S/HM/F/HM/S). The HM layers control the conversion between sin- 
glet and triplet Cooper pairing, while the F layer behaves like a spin 
filter, only allowing triplet Cooper pairs to pass across it. The F and 
HM layers are magnetically exchange coupled at the F/HM interfaces 
(B) with an angle of — n < a < n radians between the F and HM 
layer moments (a < in the figure). 



introducing an anisotropy angle a between the local magneti- 
zations of the HM and F layers. When a is or ir, the magne- 
tizations are locally parallel or antiparallel, respectively. 32 

For the strongly magnetized F layer there are two inde- 
pendent spin bands for the majority (+x) and the minority 
(— x) spins, and the singlet correlations are destroyed on a 
sub-nanometer scale. Therefore we take f° — f x = and 
identify the remaining components / ± with triplet pairs from 
the majority and minority spin bands. The governing equa- 
tions for the two independent spin bands are then 

vf cos 9 + (2lj + r- 1 ) ff = 0, (2) 

where we assume > vj . However, since scattering is 
mainly due to impurities at low temperatures, we assume that 
the pair-breaking time Tf is the same for both spin bands. 

At the interface z = z a f, between two generic layers a and b, 
the parallel component of the Fermi momentum is conserved. 
If we assume that all bands are parabolic with the same effec- 
tive electron mass, the Fermi momentum is proportional to the 
Fermi velocity, meaning f a (z ab ,6 a ,u>) = f h (z ah ,6 b ,uj) for 
v a s'm6 a = vi, sin#{,. If we further assume that v s <C Vh> v f< 
it is valid to use the single-channel approximation; the anoma- 
lous Green's function / is then only non-zero for angles close 
to or 7T in the HM and F layers, therefore cos 6 can be sub- 
stituted by ±1 in Eqs. (1) and (2). 

The phase difference between the two S layers is $ and 
the magnitude A of the bulk pairing potentials (A and Ae 1 *) 
is determined from the self-consistency equation. For T « 
Tc, we can apply rigid boundary conditions at the S/HM 
interfaces. 28 This means that / is the same at the interface as 
in the bulk of the S layers for outgoing directions. The bound- 
ary conditions at the S/HM interfaces are then f s (0, 0,oj) = 
F(A,w)l for cos 6> > and / s (d, 9, w) = F(A,w)e l * i for 
cos6> < 0, where F(A,ui) = A/VA 2 + h 2 u 2 is the equilib- 
rium value of /° and the subscript s indicates that the expres- 
sions correspond to the S side of the interface. 
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In the formalism of the LEE, the Josephson current density 
at a generic position z in the junction is given by 



J 
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/t(z,0, w )/'(z,#» 
(3) 

where 9' = ir — 9, and a z is the normal-state conductivity 
of the layer at position z. The function /' satisfies the same 
equations as / with the same boundary conditions at the in- 
terfaces, reversed exchange field {—Ik) in the HM layers, and 
exchanged spin bands (vj f-> vj) in the F layer. 



III. CRITICAL CURRENT OF THE JUNCTION 

In this section we calculate the Josephson current at the S 
side of the z — d interface and hence obtain the critical current 
of the junction. Since f s (d, 9, u) for cos 9 < is known and 
only its / component is non-zero, the Josephson current is 
determined entirely by fg(d, 9,u>) for cos 9 > 0. Due to the 
linearity of Eqs. (1) and (2), we can write 



Equations (6)— (8) were derived by setting cos6> w 1 in 
the HM and F layers. Consequently, S(lj) = S(9, lo) is in- 
dependent of 9. If we substitute Eq. (4) and f s (d,9,u) — 
, F(A,uj) e 1 * 1 (for cos 9 < 0) into Eq. (3), we recover the 
usual current-phase relation, J — Jc sin $, and the critical 
current density becomes 



f°(d,6,u) = S(9,u;)f°(0,9,uj) = S{6,w)F(A,u), 
S(6,u,) = SZ-S r S h , 



(4) 



(5) 



where the different matrices/vectors S describe the effects of 
the different layers between z — and z = d. 

The role played by the HM layers is to convert between the 
singlet component f° in the S layers and the triplet compo- 
nents / ± in the F layer, implying S/, is a 2 x 1 vector, while 
is a 1 x 2 vector. The expressions for and can be 
obtained by solving Eq. (1) in the two HM layers. The vector 
corresponding to the lower HM is 



S/, = —iBexp(—Dh) 



D 



qcosa(l — cos(5/j) + yl + q 2 sin a sin Sh 
V2(l + q 2 ) 



(6) 



(7) 



where q = Q£h, and ^ = Vh/2Ih is the correlation oscillat- 
ing length in the HM layers. We also define reduced thick- 
nesses as 8 h = d h ^l + q 2 /£ h and D a = d a (2ui + T^ v )/v a , 
where the latter one is valid for a generic layer a. Physically, 
the parameter B describes the conversion efficiency between 
singlet and triplet Cooper pairs in each HM layer, while B 2 
describes the total efficiency of both HM layers. 

Since the non-zero triplet components J* are independent 
in the F layer, S f is a 2 x 2 diagonal matrix. Integrating Eq. 
(2) gives 



S f = 



exp(-D+) 
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(8) 



therefore the 



triplet component / corresponding to the minority spin band 
decays faster than /+. 



Jc = 



dC-C ^F(A, w ) 2 Fte[SM]. (9) 

ui>0 



The integral in ( = cos 9 gives 1/2, while the sum in u re- 
quires a further approximation. Since u> <C t^ 1 for all lay- 
ers and all Matsubara frequencies with a significant contri- 
bution, we can neglect uj next to t^ 1 in D a . This implies 
that D a « d a /£ a and that S = S(oj) is independent of u> as 
well. The sum is now evaluated using J2^=o P ^+(i+2n) 2 = 
tanh {^y) and so the characteristic voltage reads 
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(11) 



The two terms in the sum are contributions from triplet pairs 
traveling through the majority and minority spin bands of the 
F layer, respectively. 



IV. DISCUSSION 

Since only triplet pairs can enter the strongly ferromagnetic 
F layer, the conversion between singlet and triplet pairs in the 
HM layers is crucial. This process is represented by the effi- 
ciency parameter B 2 in Eq. (11) and in this section we discuss 
how the conversion depends on dh, q, and a. 

We first note from Eq. (7) that B 2 is tt periodic in the 
anisotropy angle a. This implies that the critical current does 
not distinguish between locally parallel and antiparallel mag- 
netizations of the HM and F layers. In both of these cases, the 
efficiency parameter takes the form 



B 



2q 2 



'd h ^/l 



(i + g 2 ) 



2\2 



26. 



(a = 0,7r), (12) 



which shows that the periodicity of B 2 in dh is in general 



Adh = 2w^h/ v 1 + 9 2 - The periodicity Ad^ depends on the 
two intrinsic length scales of the HM layers: the oscillating 
length £/j and the wavelength A = 2ir/Q of the magnetic he- 
lix. Furthermore, Eq. (12) shows that the shorter length scale 
is the dominant one. When the oscillating length is small (i.e., 
q -C 1), we have Adh ~ 27r£/, and so the helical structure 
becomes irrelevant. The periodicity is determined by ^ as for 
uniform magnets. When the oscillating length is large (i.e., 
q ^> 1), we have Ad^ w A and so that the periodicity is de- 
termined entirely by the helical structure. It also follows from 
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FIG. 2: (Color online) The efficiency parameter B peaks at par- 
ticular HM layer thicknesses (A), and the positions and magnitudes 
of these peaks depend on the magnetic anisotropy angle a between 
neighboring HM and F layer moments (A, B). In (A) and (B), q = 1. 
For optimum HM layer thicknesses, the maximum efficiency param- 
eter -E>m ax depends sensitively on q and a (C). 



Eq. (12) that the zeros of B 2 are at dh = m Adh, and that its 
maxima are halfway between its zeros. 

If we start increasing a, the periodicity Adh remains un- 
changed, and B 2 still has zeros at dh = m Adh- However, the 
primary maxima halfway are shifted to the left, and secondary 
maxima appear with additional zeros in between (see Fig. 2). 
When the local HM and F magnetizations are orthogonal, the 
efficiency parameter reads 



B 2 = 



1 



2(1 + q 2 ) 



sin 



d h ^l 



(a = ±n/2). (13) 



The primary and secondary maxima become equivalent, the 
periodicity in dh is reduced by 2, and the zeros of B 2 occur at 
d h = mAdh/2. 

It is interesting to look at the efficiency of the singlet/triplet 
conversion at optimal HM layer thicknesses. This is deter- 
mined by q and a. In the limit when q 3> 1, the critical current 
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FIG. 3: (Color online) Experimental (bars, adapted from Ref. 19) 
and theoretical (curves) IcRn values of a Nb/Ho/Co/Ho/Nb device 
against the symmetrical thickness of Ho at 4.2 K. Theoretical curves 
are plotted for A = 3.4 nm (blue) and A = 6.8 nm (red), while the Co 
thickness is kept at 3.4 nm (note that the local magnitude of IcRn 
only weakly depends on the Co thickness). For the experimental 
IcRn values, the magnetically dead layers present at each Ho sur- 
face (w 1.2 nm) have been subtracted from the total Ho thickness. 
Theoretical curves are plotted with the following realistic parameter 



values: T c =9.1 K, v B 
ms" 1 , vj = 0.75 x 10 6 ms" 
nm, Ih = 0.25 eV, and a — 0. 



■■ 0.4 x 10 6 ms" 1 , v h = v 



= 1.0 x 10 u 



h = 4 nm, £j = 12 nm, £7=9 



vanishes. This is because if A is too small, the HM magneti- 
zation averages to zero within the length scale of In the 
opposite limit when q <C 1, we need a / for Ic ^ 0. If 
A is too large, the HM layers become uniformly magnetized, 
and we recover the results in Ref. 28 as Ic oc B 2 oc sin 2 a. 
To be more quantitative, we maximize B 2 with respect to dh', 
this calculation gives 



B, 



cos a 



sin 2 a 



2(1 



7 2 ,2 



(14) 



and the dependence of S^ ax on q and a is shown in Fig. 2(C). 
We can establish that the most efficient singlet/triplet conver- 
sion with -B 2 lax = 1/2 occurs whenever q = | cosa|. This 
is in fact an absolute theoretical maximum. Due to the two 
triplet channels / ± in the F layer, B 2 lax = 1/2 corresponds 
to perfect conversion. 

Altogether, the dependence of B 2 on the HM layer thick- 
ness shown in this paper demonstrates the profound effect that 
interfacial layer-by-layer magnetic coupling can have on the 
amplitude of the triplet supercurrent in a S/HM/F/HM/S de- 
vice. Even when the magnetizations at the HM/F interfaces 
couple parallel or antiparallel, the amplitude of the triplet su- 
percurrent oscillates with peaks and zeros commensurate on 
dh- This result agrees with experimental Ic data reported in 
Ref. 19: in Fig. 3 we directly compare the theoretical and 
experimental dependence of IcRn on Ho layer thickness and 
agreement is achieved for A in the 3.4—6.8 nm range, consis- 
tent with A values estimated in Ho thin films. 30 

This close agreement between theory and experiment pro- 
vides a clearer understanding of the role played by interfacial 
magnetic coupling between Ho and Co on the Ic found in 
Nb/Ho/Co/Ho/Nb devices. 19 The possibility of controlling Ic 
by manipulating the coupling anisotropy (i.e., a) in these or 
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similar devices has not been explored experimentally so far. 
Since the Curie temperature of Co (~ 1000 K) is significantly 
higher than either the Neel (~ 130 K) or the Curie (~ 20 K) 
temperature of Ho, it may be possible to vary a by field cool- 
ing a device from ~ 130 K. By repeating this procedure with 
the field applied at various in-plane angles, the effect of a on 
the amplitude of the triplet supercurrent Iq could be tested 
with potentially large AIc/Ic ratios obtainable. 

V. SUMMARY 

Layer-by-layer interfacial magnetic coupling between rare- 
earth helimagnets (HM) and ferromagnets (F) can strongly in- 



fluence the interconversion between singlet and triplet Cooper 
pairing in S/HM/F/HM/S-type Josephson devices. This re- 
sults in the amplitude of the spin-triplet supercurrent pass- 
ing through the F layer depending sensitively on the magnetic 
structure and the thicknesses of the helimagnetic layers. 
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